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Generation and detection of non-Abelian matrix Berry phases through manipulation 
of electric confinement potential of a semiconductor quantum dot 

S.-R. Eric Yang* and N.Y. Hwang 
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A matrix Berry phase can be generated and detected by all electric means in II- VI or III-V n-type 
semiconductor quantum dots by changing the shape of the confinement potential. This follows from 
general symmetry considerations in the presence of spin-orbit coupling terms. The resulting 2x2 
matrix Berry phase can be characterized by two numbers of geometric origin. We investigate how 
these parameters depend on the shape and area of closed adiabatic paths. We suggest how the 
matrix Berry phase may be detected in transport measurements. 
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I. INTRODUCTION 
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FIG. 1: (a) Displays electrons in quantum dots. A cyclic 
adiabatic change of the shape of the electric confinement po- 
tential does not return electron spins to the initial values, 
(b) Several adiabatic paths are shown. Adiabatic parameters 
Ai and A2 provides a means to control the shape of the dot 
electrically. 
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FIG. 2: Basic mechanism of adiabatic control is based on non- 
trivial degeneracy of II- VI and III-V semiconductors, (a) An 
electric field along the z-axis quantizes the electronic motion 
in a triangular potential along the axis. An adiabatic change 
can be induced by changing the magnitude of the electric field, 
(b) The two-dimensional electronic motion is quantized in a 
distorted parabolic potential. An adiabatic change can be 
induced by changing the magnitude of the distortion poten- 
tial. (c)In the presence of the spin-orbit terms each discrete 
eigenstate of a semiconductor quantum dot has a double de- 
generacy due to time reversal symmetry in the absence of a 
magnetic field. 



Single and few electron control by electric means in 
semiconductor quantum dots would be valuable for spin- 
tronics, quantum information, and spin qubitsp], 0, 
Adiabatic time evolution of degenerate eigenstates of a 
quantum system provides a means for controlling indi- 
vidual quantum states through the generation of non- 
Abelian matrix Berry phases 0, [y]. This method 
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can be used to perform universal quantum computa- 
tion Q. There are several semiconductor nanosys- 
tems that exhibit matrix Berry phases: they include 
excitonsQ, CdSe nanocrystalsQ, and acceptor states of 
p-type semiconductors [lOj. It is desirable to generate 
and detect matrix Berry phases by all electric means in 
semiconductor nanosystems. Recently we have demon- 
strated theoretically that it is possible to control electri- 
cally electron spins of II-V and III-V n-type semiconduc- 
tor quantum dotspj], G3 and rings [IH by changing the 
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shape of the electric confinement potential. The mech- 
anism is based on spin orbit couplings which generate 
non-Abelian vector potentials. The main ingredients of 
it are: [a] Time-reversal symmetry, which leads to dou- 
ble degeneracy [III [lij . which is depicted schematically in 
FigfT] [b] Non-Abelian U(2) gauge theory, [c] Adiabatic 
changes that break the parity symmetry of the electric 
confinement potential of the dot. When several electrons 
are present the matrix Berry phase can be generated only 
when odd number of electrons are in the dot. In such a 
case the effect of correlations can be included in a com- 
pact way (exchange effects do not affect the non-Abelian 
vector potentials) [12| . It should be emphasized that the 
presence of the matrix Berry phase follows from general 
symmetry considerations. The matrix Berry phase can 
be present even in the absence of the spin-orbit terms. 
A noteworthy example is semiconductor quantum dot 
pumps [lH, [l7j ]. which can be understood as a manifes- 
tation of a matrix Berry phase [TEl]. 

Let us give a brief explanation of matrix Berry phases. 
The electron state at time t is given by 

\y(t)) = Cl (t)\Mt)) + c2(t)\Mt)), (i) 

where |V>i.2(£)) are the instantaneous basis states satis- 
fying H(t)ipi(t) = E(t)i/)i(t) for i = 1,2 (H(t) is the 
Hamiltonian with the eigenenergy E(t) at time t). The 
matrix Berry phase $c(l) connects (ci(T), C2(T)) to 
(ci(0), 02(0)), and is given by 

<J> C (1) = ^T-pApit^dXp ^ e iY. p A p (t 1 )d\ p 

= P e i fo T, p A P dX P ^ (2) 

where time slices are ordered as t% < ... < t n < .... The 
matrix vector potentials A p 

{A P )i,i=m^), (3) 

are integrated along the path C in the parameter space 
in the order of increasing time. The path C is parame- 
terized in time as (Ai(t), \2{t)). The matrix Berry phase 
is independent of the functions Xi(t) as long as they de- 
scribe the same path C. 

In applying matrix Berry phases to II- VI and III-V 
semiconductor quantum dots there are several issues that 
need to be addressed. One issue is whether the line in- 
tegral along C can be generally converted into an areal 
integration over A in the expression for the matrix Berry 
phase. For this purpose it is useful to consider the field 
strength 2-form 

F = Y, F ^ dX ^ dX »> ( 4 ) 

{IV 

where 

Ffj,u =d„A v - d u A„ + [Ap , A v ] (5) 



with d^A v = F° r a sman area A it can be shownQ 
that the matrix Berry phase is 

$o(l) = exp(iAFi 2 ). (6) 

For an arbitrarily large A there are two cases: [1] If 
[ij,, A„] = 0, the curvature reduces to — d^A v — 
d v A^ and we can use Stokes' theorem to convert the line 
integral along C to an areal integration over A 

$c(l) = *Mi\j i^dA M A dA„). (7) 

[2] When [yl^A,] 7^ it is not possible to convert the 
line integral into an areal integration. We will investigate 
the applicability of Stokes' theorem in the case of Rashba 
and Dresselhaus terms. 

Second issue is that it is often complicated to calculate 
the non-Abelian vector potentials, and, consequently, it is 
difficult to compare experimental and theoretical results. 
This is especially so when many electrons are present. 
It would be useful to have a simple, yet general, mathe- 
matical expression for the matrix Berry phase. We will 
show that the matrix Berry phase is a unitary 2x2 ma- 
trix and that it may be thought of as a spin rotation 
matrix |18|. A spin rotation matrix is characterized by 
the direction and angle of rotation. It is unclear how the 
geometric properties of the adiabatic path are reflected 
on them. We show in this paper how the information of 
the path C can be encoded in the direction and angle of 
the rotation. 

Another issue is whether the matrix Berry phase can 
be detected in a transport measurement, which would 
complement the detection of matrix Berry phases by in- 
frared optical measurements (ll) . 
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FIG. 3: The geometric information about the closed curve 
in the parameter space is encoded into the rotation axis rh 
and angle 2a. When only the Rashba term is present m is 
independent of path C. 



Let us now give a brief summary of the main results of 
this paper. First, the matrix Berry phase has the form 

$ c (n) = e ii^n)rn,-S = cos ( an ) J + j sin ( an ) ? 7i . ^ ( 8 ) 

where a are Pauli spin matrices and n is the number 
of periodic adiabatic cycles. The information about the 
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closed curve in the parameter space is encoded into the 
rotation axis 



m = (Re(/3), -Im(/3), y/T^W) (9) 

and the angle of rotation 2cm. The constants a and (3 
are real and complex numbers, respectively. This result 
is even valid for many electron systems with correlations 
as long as double degeneracy is present. The expression 
Eq.® is valid only at time t = nT, and it should be 
stressed that in the time interval (n — 1)T < t < nT 
the probability amplitudes c\(t) and c^it) exhibit a much 
more complicated behavior, see Figj7l Second, we find 
(3=1 and in = (1,0,0) when only the Rashba term is 
present, and all the geometric information is contained 
in the rotation angle a. We have derived an analytical 
expression for a. Furthermore, remarkably [A^, A u \ = 
and Stokes' theorem can be applied. When both Rashba 
and Dresselhaus terms are present each path has different 
m and a. We have investigated how they depend on the 
shape of each path. Moreover, we find that L4 M , A v \ ^ 
and that Stokes' theorem cannot be applied. Third, we 
propose how the presence of the matrix Berry phase can 
be measured in transport experiments by changing the 
electric confinement potential suddenly, see Sec. III. This 
is shown schematically in FigdJ 
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FIG. 4: (a) A matrix Berry phase is generated when the 
adiabatic parameters Ai and A2 go through a cyclic change. 
(b)After a matrix Berry phase is generated the Rashba elec- 
tric field is suddenly increased so that an electron can tunnel 
out the dot. 



II. SEMICONDUCTOR QUANTUM DOTS 
WITH THE RASHBA AND DRESSELHAUS SPIN 
ORBIT COUPLINGS 

A. Model 

The Hamiltonian is H = + U{x,y) + V(z) + 

+ Hj). We take the two-dimensional potential to be 
U(x, y) = ^m*oj 2 x 2 + ^m*u> 2 y 2 + e'y, where the term e'y 
represents a distortion of the two-dimensional harmonic 
potential. The strengths of the harmonic potentials are 
denoted by lu x and io y . The characteristic lengths scales 

along x- and y-axis are R x , y = \J ~- hi our work the 

triangular potential V{z) is sufficiently strong and only 
the lowest energy subband is included. The characteris- 
tic length scale along the z-axis is 1/ y // 0.8(2meE/h 2 ) 2 ^ 3 , 
where E is the Rashba electric field applied along the 
z-axis. The Rashba spin orbit term[19|, |2fJ is 

Hr = cr (o- x k y - o-yk x ) , (10) 

and the Dresselhaus spin orbit term[2l[ is 

H D = c D ((a x k x (k 2 y - k 2 z )) + (a y k y (k 2 z - k 2 x ))) . (11) 

Note that the Rashba spin orbit constant cr depends on 
the external electric field E. Here k x ,y are momentum 
operators (k x = and similarly with k y ). The Hamil- 
tonian matrix is represented in the basis states \mno~) of 
the harmonic oscillator states in the xy plane with spin 
component a. The subband wavefunction f(z) is sup- 
pressed in the notation \mna). Let us take uj x = 2u> y 
(Other values of u x can also be chosen). Then the low- 
est eigenenergy state of H is |mn) = |00) with the en- 
ergy Eq = %hu) y and the next lowest eigenenergy state 
is 1 01} with the energy E\ = |i?o- The typical value 
of the energy spacing between the quantum dot levels, 
Eo, is of order 1 — lOmel^. It can be several times 
larger in self-assembled dots. The energy scale of the 
Rashba term is Er = cr/R ~ 0.01 — lOmeV, where 
the length scale R ~ 100A is the lateral dimension of 
the quantum dot. The energy scale of the distortion 
potential is E p = (Q\e'y\i). Its magnitude is of or- 
der 1 — lOmeV, depending on the electric field applied 
along the y-axis. The energy scale of the Dresselhaus 
term is Ejy = cp/R 3 , and it can be larger or smaller 
than the Rashba term, depending on the material (l9j. 
Here E D = ic D {Q\k y \l){{f{z)\k 2 z \f{z)) ~ (0|Ag|0» = 
cd/R v R 2 — cd/R v R 2 , and it originates from the second 
term of the Dresselhaus term (the first term in the Dres- 
selhaus term is zero in our model since (0|fca;|0) = 0). 

We will use a truncated version of the Hamiltonian 
matrix, which makes it possible to write the eigenstates 
as a linear combination of four basis states made out of 
\mn) and spin degree of freedom: 

W) = c OAT |0,0,T>+c o ,i, T |0,l,T} 

+ c 0) o,i|0,04)+co, u |0,U). (12) 
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The advantage of this truncated model is that it is exactly 
solvable. In this paper we choose the following eigenstates 
of the lowest energy shell as the instantaneous basis vec- 
tors: 



3E V 
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(13) 



The non-Abelian vector potentials are evaluated with re- 
spect to these basis vectors. 

In the following we choose the adiabatic parameters as 
the Rashba energy and the distortion energy 



Ai = E R = c R /(V2R y ) 
A2 = E„ — 



(14) 



The cyclic adiabatic path is given by E R (t) = E RjC 
AE R cos(ujt) and E p (t) = E p<c + AE p sm(ujt). 



B. Geometric dependence 



When only the Rashba term is present the structure 
of the non-Abelian vector potentials is remarkably sim- 
ple. The orthonormalization = &y gives that the 
diagonal matrix elements (A p )i t i are real and that the 
off-diagonal elements satisfy {A p )ij — (Ap)*^. We find 
non-Abelian vector potentials A\ and A2 
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A, = 
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El) 



(15) 



(16) 



(Note that (f(z)\g§-\f(z)) = 0.) The diagonal elements 
An are zero and the off-diagonal elements Aij are real. 
From this it follows that the parameter (3=1 and only 
one geometric number a is needed. The matrix Berry 
phase is given by 



$c(«) 



_ / cos(cm) i sin(cm) 
i sin(an) cos(cm) 



(17) 



One can show that [Ai,^] = when only the Rashba 
term is present, and that Stokes' theorem can be applied, 
Eq. ([7]). The field strength has a simple form F 12 = 
fa x jEl where 



AE„ 




FIG. 5: For elliptic paths with fixed area with AEr changing 
from 3.2Eo to 4.8i?o. In the presence of both of the Dressel- 
haus and Rashba terms (solid line). Only the Rashba term 
present (dashed). How j3 changes with semiaxis AEr in the 
presence of both of the Dresselhaus and Rashba terms. 



/ = - 



2 (9x 2 + 9y 2 + 1) 



3/2 ' 



(18) 



with x = E p /E and y = E R /E . From this we can 
compute the rotation angle 



fdxdy, 



(19) 



which gives \a\ < 0.7854. 

When both the Rashba and Dresselhaus terms are 
present Stokes' theorem is not applicable. We have inves- 
tigated how m and a depend on the shape of the paths. 
They can be computed numerically by solving the time 
dependent Schrodinger equation 



ihdi 



3 



(20) 



The matrix elements A^ are given by Ai 



where the sum over p in A^ is meant 



to be the sum over X p . We have considered elliptic paths 
with the area ttAE r AE p = 16-kEq with E RtC = 5E and 
E PtC = 5Eq. The results are shown in FigfS] (Note that 
the dependence of Ed on the adiabatic parameter E R is 
not well-known, but the essential physics does not depend 
on the exact functional form, as discussed in Ref . [Ill] ) . 
FigJS] also shows numerical results for a when only the 
Rashba term is present. We have verified that these nu- 
merical results agree well with the results obtained from 
the equation Eq. (fT9"|) . 
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C. Time dependence 




(EM t 



FIG. 6: Probabilities as a function of time in the presence of 
the Rashba spin orbit coupling. Black dots are fitted values 
at t — nT using the analytical expression for the matrix Berry 
phase, Eq.JTTJ). 



Let us first investigate the time evolution of the system 
in the presence of the Rashba term. The time-dependent 
Schrodinger equation is solved using the Runge-Kutta 
method. For the parameters Eb,, c = 2£ , o 5 Ep,c = Eq, 
AE r = 1.9E , AE p = 0.9E , and frw Y = O.LE we 
calculate |ci(t)| 2 and |c 2 (i)| 2 , see FigJHl We can fit the 
data to \c\{nT)\ 2 = cos 2 (an), \c2(nT)\ 2 = sin 2 (an), with 
a w 0.182. In this case the matrix Berry phase takes a 
simple form with (3—1. We have tested our numeri- 
cal method as follows: For two different periods T-y and 
T 2 we find the matrix Berry phases are the same, i.e., 
c 2 (Ti) = c 2 (T 2 ) and ci(T x ) = Ci(T 2 ). This is consistent 
with the fact that the matrix Berry phase is a geomet- 
ric effect independent of how the path is parameterized. 
Let us now include the Dresselhaus spin orbit term as de- 
scribed in Ref . In this case the matrix Berry phase is 
more complicated since (3^1. The diagonal elements An 
are non-zero and the off-diagonal elements A,j are com- 
plex numbers. Using the same parameters -Er jC , E PiC , 
AEr, AE p , and u> as in the previous calculation, we 
have computed | c x (t) | 2 and | C2 (i) | 2 , which are displayed 
in Fig|7] These results can be fitted with the parameters 
a w -0.477 and (3 « -0.190 + z0.965. In both cases the 
numerical results are consistent with those of the analyt- 
ical expression at t — nT, given by Eq.®. In the time 
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FIG. 7: Probability c\(i) when Dresselhaus spin orbit inter- 
action is included. Black dots are fitted values using the an- 
alytical expression for the matrix Berry phase, Eq.([8j|. The 
initial state is Ci(0) = 1 and £2(0) = 0. 



interval (n — 1)T < t < nT the coefficients ci display 
a complicated behavior. It should be stressed that even 
when the adiabatic path is not closed the actual state 
^(t) is given by a nontrivial linear combination of the 
instantaneous eigenstates ipi(t) and "02 (t). In this case a 
matrix phase of the form 

$ c = Pe*^S I ,V^ i (21) 
will appear, where the path C is unclosed. 
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III. DISCUSSIONS AND CONCLUSIONS 

We have investigated II- VI and III-V n-type semicon- 
ductor dots with spin orbit coupling terms and have 
shown in this paper that the geometric information of a 
closed adiabatic path can be encoded in the rotation axis 
and angle of the spin half matrix. This result is appli- 
cable even for strongly correlated states, as along as the 
states are doubly degenerate. In our truncated model for 
the Hamiltonian the matrix Berry phase can be found 
exactly when the Rashba term is present. In this case 
the Stokes' theorem is applicable and the relevant field 
strength can be calculated analytically. 

Even when the adiabatic path is not closed the ac- 
tual state,^(t), at time t will not be the instantaneous 
eigenstate ipi (t) . Instead it is a linear combination of the 
instantaneous eigenstates 

|*(t)) = ci(t)|^i(t)) + ca(t)hk,(t)). (22) 

The expansion coefficients Ct(t) change in a non-trivial 
way in the interval (n — 1)T < t < nT , see Figs|6]and[7] 
Experimental investigations of matrix Berry phases in 
II- VI and III-V semiconductor quantum dots would be 
most interesting. Optical dipole measurements [22j were 
proposed before 11]. Here we propose a transport mea- 
surement. In order to measure the matrix Berry phase 
we perform the following set of manipulations: 

1. Prepare an initial state by applying a small mag- 
netic field and taking the zero field limit, as ex- 
plained in detail in Ref. [TT| | . 

2. Perform an adiabatic cycle by changing En and E p 
that defines the shape of the electric confinement 
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3. Change suddenly the electric field along the z-axis 
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FigHb). 
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current Q. 

5. Repeat this set of procedures many times. 
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